Sussex preprint SUSSEX-AST 95/12-1, 



istro-ph/9512102 



Cold dark matter models with a cosmological constant 



Andrew R. Liddle, 1 David H. Lyth, 2 Pedro T. P. Viana 1 and Martin White 3 

1 Astronomy Centre, University of Sussex, Falmer, Brighton BN1 9QH 

2 School of Physics and Materials, University of Lancaster, Lancaster LAI J^YB 

^Enrico Fermi Institute, University of Chicago, 5640 S. Ellis Ave, Chicago, Illinois 60637, USA 



3 



> 

O 

(N 

in 

Oh 
i 

o 



Accepted 1996 April 29. Received 1996 April 19; in original form December 18 



ABSTRACT 

We use linear and quasi- linear perturbation theory to analyse cold dark matter models 
of structure formation in spatially flat models with a cosmological constant. Both a 
tilted spectrum of density perturbations and a significant gravitational wave contri- 
bution to the microwave anisotropy are allowed as possibilities. We provide normal- 
izations of the models to microwave anisotropics, as given by the four-year COBE 
observations, and show how all the normalization information for such models, in- 
cluding tilt, can be condensed into a single fitting function which is independent of 
the value of the Hubble parameter. We then discuss a wide variety of other types of 
observations. We find that a very wide parameter space is available for these models, 
provided Qq is greater than about 0.3, and that large-scale structure observations show 
no preference for any particular value of f2o in the range 0.3 to 1. 

Key words: cosmology: theory - dark matter. 



1 INTRODUCTION 

Even before the COBE satellite's detection of microwave 
background anisotropies, a popular strategy for dealing with 
the perceived ills of the standard cold dark matter (CDM) 
model of structure formation was to reduce the matter den- 
sity. This shifts the redshift of matter-radiation equality, dis- 
torting the shape of the spectrum of density perturbations to 



provide a good fit to the galaxy correlation function. In order 
to retain compatibility with typical models of inflation, a fiat 
universe is desired which can be brought about by adding a 
cosmological constant A to the model (Peebles 1984; Turner, 
Steigman & Krauss 1984; Efstathiou, Sutherland & Maddox 
1990; Kofman, Gnedin & Bahcall 1993; Stompor, Gorski & 
Banday 1995; Klypin, Primack & Holtzman 1995). We will 
denote such models ACDM models. 

The ACDM models are most commonly considered in 



mordial spectrum of density perturbations ( 


Efstathiou et al. 


199C; 


Stompor et al. 1995 


; Klypin et al. 1995 


), corresponding 



to spectral index n — 1. Another limitation usually placed 
on the models is that the gravitational waves should not 
contribute significantly to the large-scale microwave back- 
ground anisotropy. However, within a given model of infla- 
tion both the spectral index and the gravitational wave con- 
tribution are determined, and while many models do give 
n very close to 1 and negligible gravitational waves, there 
are other models that do not. It is therefore realistic to allow 
both n and the gravitational wave contribution to vary when 
confronting a model with the data, and we will adopt that 



viewpoint in this paper. Observational constraints on n and 
the gravitational wave contribution will in the future place 
powerful constraints on models of inflation, and hence on the 
nature of the fundam ental interactions at very high-energy 
scales. Kofman et al. (1993) have previously considered one 



example of the model with (mild) tilt and a gravitational 
wave contribution. Because the ACDM model has generated 



further interest r ecently (Krauss & Turner 1995; lOstriker & 
Steinhardt 199£ ; Bagla, Padma.iiabhan fc Naiiikar 1995| ), wc 
believe that the full parameter space merits detailed study, 
and in this paper we compare a variety of predictions of the 
model, accessible via linear and quasi-linear perturbation 
theory, with observations of various types. 

We thus permit arbitrary variation of three main pa- 
rameters, namely the matter density Q.q expressed as a frac- 
tion of the critical density (the assumption always being that 
a cosmological constant restores spatial flatness) , the present 
value of the Hubble parameter Ho = 100 h kms~ Mpc~ , 
and the spectral index n of the primordial spectrum. In ad- 
dition, we will consider two possibilities concerning gravi- 
tational waves; first that they are negligible, and secondly 
that they have an amplitude appropriate to the power-law 
inflation model. 

One of the most important observations is that of 
cosmic microwave background (CMB) anisotropies by the 
COBE satellite, and we provide careful normalizations of 
the models to the four-year COBE data. We quote a new 
fitting function which gives accurate normalizations for any 
combination of parameters. We also provide some compar- 
ison with intermediate-scale anisotropies, where the obser- 
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vational situation is less clear. We then go on to address a 
variety of constraints concerning large-scale structure, ex- 
tending techniq ues we have already employed for critical- 
density models ( Liddle et al. 1996b) and for open universe 
CDM models ( Liddle et al. 1996a ). The most closely related 
paper is Stompor et al. ( 1995[ ), who provided extensive com- 
parison with COBE for the case of a scale-invariant (n = 1) 
initial power spectrum, along with discussion of other large- 
scale structure constraints. They did not, however, discuss 
tilt or gravitational waves; nor did they attempt to delineate 
the allowed parameter region in the flo~h plane. 



2 THE POWER SPECTRUM 

We specify the power spectrum of the densit y contrast, fol- 
lowi ng the notation of Liddle & Lyth ( 1993| ) and Liddle et 
al. ( |L996a| ), as 



*M*) = (^) S Uk)T 2 (k) 



!, U rf, u ff 2 (») 



(i) 



where k is the wave number, a is the scale factor, H = a/a is 
the Hubble parameter and Q is the density parameter. Sub- 
script '0' indicates the present value. One can think of the 
remaining terms in three parts. First, S 2 i(k) gives the pri- 
mordial spectrum of anisotropies, before they are affected by 
any evolution. The shape is determined by whichever mech- 
anism is used to create the spectrum, which we will assume 
is cosmological inflation. If S 2 i (k) is a constant then one has 
the Harrison-Zel'dovich spectrum, while 5 2 I (k) oc fc n_1 cor- 
responds to a 'tilted' spectrum with spectral index n. S uch a 



spectrum is the g eneric prediction of slow-roll inflation ( Lid- 
dle & Lyth 1992). The amplitude of S 2 l (k) fixes the present- 



day normalization of the spectrum. The transfer function 
T(k) contains the evolution of the spectrum from its primor- 
dial form to its present form, and depends on cosmological 
parameters and the nature of the matter in the universe. 
It approaches unity on large scales (which retain their pri- 
mordial form), and in a CDM model is time-independent 
well after matter-radiation equality. Finally, there is a fac- 
tor which specifies the time, or redshift, dependence of the 
amplitude of the spectrum at late epochs. If Qq equals 1 
this is carried in the {aH) 4 prefactor; for low-density cos- 
mologies there is an additional suppression factor which we 
parametrize via a quantity g(0,), defined below. 

A quantity related to the power spectrum is the dis- 
persion <j(R) of the density field smoothed on a (comoving) 
scale R, defined by 



a 2 (R) 



W 2 (kR)V s {k) 



dfc 

T 



(2) 



To carry out the smoothing, we will always use a top-hat 
window function W(kR) defined by 



W(kR) = 3 



s'm(kR) cos(fc.R)\ 
(kRf {kR) 2 ) 



(3) 



For CDM models as considered here, t he tr ansfer func- 
tion is accurately given by Bardeen et al. (1986) as 



7cDM(g) 
[1 



In (1 + 2.34a) 

— - x 

2.34g 

+ 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71q) 4 l ~ 1/4 



(4) 



with q — k/hT, w here the 'shape parameter' T is defined as 
( [gugiyama 1995 ) 

r = n h exp(-n B - n B /n ) . (5) 

Here $1b is the baryon density, which we take to equal 
O.OI6/1 -2 as suggested by recent analyses of nucleosynthe- 
sis (Copi, Schramm & Turner 1995a, b). This fit to V is good 
for both flat and open universes (unless Qq is very low). 

We have accounted for the redshift dependence of the 
amplitude of the power spectrum by introducing a growth 
s uppr ession factor g(Q), following Carroll, Press & Turner 
( [l992| ). This gives the total suppression of growth of the 
dispersion u(R) relative to that of a critical-density universe, 
and is accurately parametrized by 



1 209ft 
70 + 140 



140 



+ D. 



4/7 



(6) 



This formula can be applied at any value of ft. For a matter- 
dominated flat universe, the redshift dependence of fl is 
given by 

A3 



q(z) = n 



1- 



i - n + (i + z) 3 n 



(7) 



Since the growth law in a critical-density universe is a(z) oc 
(1 + z)' 1 (carried by the (aH) 4 term in equation ([j])), the 
redshift dependence for arbitrary Qq is therefore 



a(R,z) = a(R,0) 



g("(*)) 

9 (n ) 



(8) 



3 MICROWAVE BACKGROUND 
ANISOTROPIES 

3.1 COBE 

With the detection of large-angle temperature fluctuations 
in the cosmic microwave background, the COBE satellite 
made possible for the first time an accurate normalization of 
models of structure formation. At the large scales probed the 
spectrum can be normalized in a regime where the theory is 
well understood, circumventing complications introduced by 
processing of the primordial spectrum and the relationship 
between observed structure and the underlying density field. 

The first-year COBE data were of low signal-to-noise ra- 
tio, with the rms fluctuation having a 30 per cent error. Fits 
to the full data set, the angular correlation function or the 
rms fluctuation all gave consistent va lues for the quadrupole 
normalization: Q rm s-ps = 17=1= 5 /iK (Smoot et al. 199S ; |Scl 



et al. 



1994p . 



ak & Bertschinger 1993; ; |Scaramella fc Vittorio 1993j ; [WiiglitJ 



The second year of data (Bennett et al. 1994) resulted 
in a dramatic improvement of the signal-to-noise ratio and a 
consequent increase in the degree of refinement of the anal- 
yses. The two-year data constrain the large-scale normal- 
ization to within 10 per cent (at la), with 5-7 per cent of 
this being due to irremovable cosmic and sample variance. 
Ironically, along with the better signal-to-noise ratio came 
ambiguity in the number to use for normalization (amount- 
ing to a 30 per cent discrepancy!), due in large part to a 
low quadrupole in the two-year map (Banday et al. 1994; 
Bunn, Scott & White 1995). Normalization of models di- 
rectly to the temperature maps became essential to obtain 
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all the in formation a v ailable from the COBE data ( Gorski et 
al. 199^; |Bond 1995] ; |Bunn 1995[ ). In addition, highly accu- 
rate calculations of theoretical predictions and their relation 



1995a |; Korski et al. 1995|~ [Tegmark fc Bunn 1995|; [White & 



to large-scale structure were include d in the analyses ( Bunn 
et al. 1J)95{ [Bunn & Sugiyama 1995|; |Hu, Bunn & Sugiyama 



Bunn 1 | 995 ; Btompor et al. 1995; Yamamoto fc Bunn 1995^ 



Cayon et al. 1996; White & Scott 1996), making the COBE 



normalization one of the most accurately known pieces of 
information about large-scale structure. 



W ith the release of the four- year anisotropy maps ( Ben- 
nett et k\. 1996;: |Banday et al. 1996] ; |G6rski et al. 1996| ; |Hm 
shaw et al. 1996), we have in hand all of the knowledge about 



large angular scale anisotropics that we can expect to ob- 
tain in the near future. The full data set prefers a slightly 
lower normalization (~ la) than the two-year data, due in 
equal parts to a stat istical downward flu ctuation and an im- 
proved galaxy cut (Gorski et al. 199C). Models with near 



scale-invariant spectra, such as the ACDM models consid- 
ered here, are now a better fit to the data and the actual 
quadrupole on the sky is no longer anomalously low. How- 
ever, while the 'final' normalization is lower than that of 
the two-year data, the central value is still higher than that 
which would be obtained from considering only the rms fluc- 
tuation in the map. The COBE data still cannot be summa- 
rized by one number and fits to the full data set are necessary 
to obtain a precise normalization. 

For low-density CDM models, the COBE data play a 
vital role in breaking the degeneracy between flat ACDM 
models and open CDM models. While arguments based on 
structure formation alone are insens i tive to a cosmological 
constant ( |Martel fc Wasserman 199C ; Martel 1991; Lahav et 
199l|), the COBE normalizations of the two theories are 



al 



very different. We will say more on this in the conclusions. 

Inflation predicts not only a spectrum of density per- 
turbations, but also one of gravitational waves (also known 
as tensor perturbations). Those with a wavelength compa- 
rable to the size of the observable Universe can induce extra 
microwave background anisotropics over and above those 
caused by the density perturbations. In absolute general- 
ity, the contribution of these to the COBE normalization is 
independent of the spectral index of density perturbations, 
permitting an arbitrary reduction in the normalization of 
the density perturbations. For n < 1, we choose to examine 
two cases. The first case has negligible gravitational waves, 
which is a prediction of many inflationary models even when 
they give significant tilt (the reason being that the gravita- 
tional waves are negligible unless the potential during in- 
flation is within a few orders of magnitude of the Planck 
scale). In the second, the gravitational wave amplitude is 
the one given by power-law inflation (PLI); this happens to 
be a good approximation for the known models of inflation 
that give a significant gravitational wave contribution. We 
will also study n > 1, but in that case we will not consider 
gravitational waves since it is hard to arrange gravitational 
waves of significant amplitude when n is greater than 1. 

For our assumed inflationary theories the density per- 
turbation and gravitational wave spectra are given by per- 
fect power laws. We compute the corresponding radiation 
power spectrum by numerically integrating the coupled Ein- 
stein, fluid and Bolt zmann equations to the present as dis- 
cussed in Hu et al. (1995b). We carry out a fit to the four- 



year COBE data using the method of White & Bunn (199 
with the customized galactic cut of Kogut et al. (1996). A 



full description of this calculation will be given by Bunn & 
White (in preparation). 

Defining the spectrum in the manner of equation ([]]) is 
extremely useful, because when expressed in terms of <5h the 
normalization is independent of the present Hubble param- 
eter h to an excellent approximation. Further, the scale of 
the COBE observation is large enough that details of the 
material content of the Universe are unimportant, making 
the normalization independent of the baryon density and 
whether or not the dark matter is multi-component. So the 
present-day normalization of the power spectrum depends 
only on the density parameter Qq, the tilt of the spectrum 
n and, if included, the amplitude of any gravitational waves 
that might influence the COBE result. 

In the scale- invariant case 6n(k) is constant, and can 
be specified on any scale. More generally, it is necessary to 
specify the scale on which it is quoted, via 



5l{k) = 8l{k Q ) (|-) 



(9) 



We will choose ko = a^Ho, the present Hubble scale, and 
write 8h = fo(fco). 

We find that we are able to represent the normalization 
of tilted ACDM models by a single fitting function 



&(n,n ) = 1.94 x 10" 5 ft" 



0.785-0.05 ln!in 



where 
f(n) 



-0.95(n 



0.169(n 



= 1.00(n - 1) + 1.97(n - 1)' 



exp [/(n)] , (10) 



No tensors ; (11) 
PLI. (12) 



These fits are accurate to 3 per cent in the range 0.2 < 
Qo < 1 and 0.70 < n < 1.2 (the latter formula though only 
applying for n < 1). The statistical uncertainty in all of 
these numbers is 7 per cent, and there is an additional 3 per 
cent 'systematic' uncertainty from the choice of galaxy cut; 
we combine these to obtain a 15 per cent uncertainty at 2a. 
The fits are accurate for any reasonable h and J1b, and for 
any choice of dark matter content. This extremely simple 
form is much m ore ea sily applied than th e tabu lar data of 
White & Bunn (1995) and Stompor et al. (1995), as well as 



applying to a much wider parameter space of models. 

Although we normalize all models to the central COBE 
value, we allow for its uncertainty by adding it in quadrature 
to the relative error of other observations constraining the 
amplitude. 



3.2 Intermediate-scale anisotropics 

Potentially the most powerful test structure formation mod- 
els will have to surpass in the future is the shape of the 
microwave anisotropy power spectrum on smaller angular 
scales than those sampled by COBE. It is expected that 
this can be measured to great accuracy by proposed satel- 
lite observatories. At present the observational situation is 
rather uncertain, and so it is best to focus on a single quan- 
tity — the height of the peak at degree scales relative to the 
large angular scale plateau. 

In a ACDM model, the first peak in the spectrum is 
located at around £ ~ 220, with a small dependence on Oq 
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and h. Here £ is the multipole number in an expansion of 
the CMB temperature field: AT = ^2 tm atmYe m (9 , <f)) with 
Yi m (8,(f)) the usual spherical harmonics. The COBE data 
sample multipoles I = 2-30, with I ~ 10 being the pivot 
point with changing n. We define a quantity 



£(£+l)Ce 
HOCio 



(13) 



where the Ce are (as usual) the expected values of the coef- 
ficients of the multipole expansion: Ci = {\on m \ 2 ^. Focusing 
on D(220), th e predicted v alue is given by the following fit- 
ting function ( [White 1996| ): 



0(220)=; 5.1 (^2)", 



(14) 



where we have formed the degenerate combination of pa- 
rameters 

v = (n - 1)- 0.32 ln(l + 0.76r) +6.8(fi B h 2 - 0.016) 

-0.371n(2/i) - 0.161n(fi ) - 0.65r . (15) 

Here r is the tensor to scalar ratio, normalized to r — 7(1— n) 
for power-la w infl ation in the fio — 1 and n — > 1 limit as in 
Davis et al. (1992). This form accurately describes the peak 
height for Q > 0.3 and 0.01 < n B h 2 < 0.02. For models 
with lower f2o the shape of the CMB spectrum is sufficiently 
different from when f2o = 1 that a case by case comparison 
is warranted. 

The term depending on r shows the reduction in the 
peak height due to reionization between z = and the last 
scattering surface. The optical depth to Thomson scattering 
for ionized fraction x e from z = to the reionization redshift 
zr is 



\M)(l + ZR.) 3 + l-fi0 



(16) 



r = 0.035 ^ hie 

"0 

The uncertainties in the epoch of reionization and value of 
fish 2 are the biggest barrier to using the measured height 
of the peak to constrain the spectral slope. 

Present observations, focusing on the COBE data and 
intermediate-scale experiments which probe only around the 
peakQ, provide the constraint v = — 0.02±0.12 (±2<j range). 
While this result provides a strong constraint in the case of 
the £lo = 1 models, limiting the degree to which tilt can be 
used to reduce small-scale power, it is not very constraining 
for ACDM models. The reason for this is that ACDM mod- 
els have a larger amount of power on degree scales than the 
critical models, giving them more room for a tilt. The up- 
per limit on the height is compromised by the possibility of 
reionization in these models, which reduces the degree-scale 
power and is exponentially sensitive to the assumed redshift 
of reionization. We will only use the lower limit, and make 
the conservative assumption of no reionization. 

We have fixed the baryon density in equation ( p^| ) to the 
central nucleosynthesis value. We include the uncertainty 



We include the new Saskatoon-95 data (Ncttcrficld et al. 1996) 



in this analysis, marginalizing over the 14 per cent calibration 
uncertainty. We also include the COBE, Python III 3-point, SP94, 
MAX, and MSAM three-point da ta wh ich probe predominantly 
around the first peak; see White (1996) for a discussion of these 
experiments. We do not perform any foreground subtraction in 
this analysis. 



in fish 2 by adding it in quadrature to the observational 
limits on v, though this only increases the error bar on v to 
±0.13. This uncertainty would become greater, however, if 
one were to allow Q B h 2 to drift outside the nucleosynthesis 
range. We note in passing that the next generation of smaller 
angular scale (mostly interferometer) CMB experiments will 
allow tighter constraints, and more importantly will suffer 
less from uncertainties in cosmological modelling. 



4 LARGE-SCALE STRUCTURE 
OBSERVATIONS 

Our use of observational data has been described exten- 
sively for the case of a critical-density universe in Liddle 
et al. (19961), and we will not reproduce the analyses here 



except where important differences arise when the cosmo- 
logical constant is introduced. Instead we provide a short 
summary. 

4.1 Galaxy correlations 

We use results fro m the compilation of data carried out by 
Peacock & Dodds (1994). They supply values of the power 
spectrum at a selection of different k, which can be used 
directly to constrain the shape parameter P. In their paper 
they quoted a central value of r = 0.255 + 0.32(l/n — 1 ); 



however, the y have noted since ( Peacock 199t ; Peacock & 
Dodds 1996) that both non- linear bias and non- linear evo- 
lution of the power spectrum may have some effect on short 
scales. We choose to drop the four shortest scale points and 
refit. We find P = 0.230 + 0.28(l/n - 1), where the un- 
certainty is plus 18 per cent and minus 15 per cent at 95 
per cent confidence. The galaxy correlation function data 
are not useful for constraining the amplitude of the power 
spectrum. 

4.2 POTENT 

The bulk velocity flow has long been suggested as one of the 
main worries for the low-density model, which predicts lower 
velocities than its critical-density counterpart for the same 
size of density contrast. However, the COBE normalization 
for low densities is actually much higher, as we have seen, 
and it will turn out that the velocities give no useful con- 
straint. In fact, at present the bulk velocity constraint is the 
weakest constraint of all. Additional ways of using the ve- 
locity data, which permit stronger constraints not directly 
interpretable in linear theory, will be discussed in Section 

We use the measurement of the bulk flow around our 



posit ion by the POTENT technique ( Bertschinger fc Deke] 



198E), using the Mark III data set (Dekel 1994). This con- 



tains the evaluation of the bulk flow in spheres about our 
position for a variety of radii. However, the bulk flow is sen- 
sitive to long wavelengths to a much greater extent than 
the dispersion of the density field, and consequently in any 
given realization these data are highly correlated. Because 
of this, we concentrate on the measurement on a single scale, 
40 ft -1 Mpc. This is obtained observationally first by carry- 
ing out a smoothing with a 12 h^ 1 Mpc Gaussian to generate 
a continuous velocity field, and then velocity reconstruction 



© 0000 RAS, MNRAS 000, 000-000 



CDM with a cosmological constant 5 



is carried out to give the flow corresponding to a top-hat 
smoothing. The theoretical prediction for the rms bulk flow 
is therefore 



V(fci?)e-< 12h ~ lfc > 2 ^, (17) 



where W(kR) is the top-hat window given by equation 
and the factor / 2 (f2o) gives the velocity suppression. Often 
/(fio) is approximated as f2o' 6 , which turns out to work well 
for both open and flat models (Lahav et al. 1991), but for 



the sake of accuracy wc numer ically obtain the precise result, 
following Carroll et al. (|l992[). 



The Mark III POTE NT analysi s gives the bulk flow in 
a 40/i _1 Mpc sphere as (|Dekel 199^) 



«obs(40 /i _1 Mpc) = 373 ± 50 kms" 



(18) 



The quoted error represents different strategies for coping 
with sampling gradient bias and is to be considered as an 
estimate of the systematic error. In addition there is a la 
random error of 15 per cent. However, both of these are dom- 
inated by cosmic variance, which arises because we have only 
a single measurement which is to be drawn from an ensem- 
ble following a X3 distribution. Modelling this al ong with 
the observational errors as in Liddle et al. (1996b), we find 



that the 95 per cent confidence limits on the estimate of the 
dispersion are +295 per cent and —47 per cent. The huge 
asymmetry is caused by the asymmetry of the xi distribu- 
tion, and only the lower limit is of interest to us. 



4.3 Cluster abundance 



In a recent paper, Viana & Liddle ( 1996 ) re-analysed the 
constraint on the power spectrum from the abundance of 
large galaxy clusters, using X-ray observations and a calcu- 
lational technique based on Press-Schechter (1974) theory. 
Traditionally, the cluster abundance constraint is quoted on 
the scale 8 ft -1 Mpc, which in linear theory corresponds more 
or less to the appropriate mass of a cluster. The result found 
was 



a s = 0.60 fi 



C(fJ ) 



where the fitting function 
C(fio) = 0.59 - 0.16 fi + 0.06 ft 2 , 



(19) 



(20) 



parametrizes the changing power-law index of the fto de- 
pendence. At fio = 1, the uncertainty was found to be +32 
per cent and —24 per cent. As flo is decreased, the uncer- 
tainty becomes slightly larger; this increase can be expressed 
by multiplying the uncertainties for fio = 1 by a factor 
^o.26io glQ o ^jjQygjj difference does not really become 
important until Qq < 0.3). 



4.4 Early object formation 

For critical-density models of structure formation, such as 
models with both cold and hot dark matter, an important 
constraint is whether or not structure can form sufficiently 
early. The strongest constraints available at the moment 
come from the amount of gas in damped Lyman alpha sys- 
tems, which for those models imposes new constraints on 
parameter regions not excluded by other data. 



For low-density models, whether open or spatially flat, 
these constraints are much less of a concern, for the rea- 
son that structure has been growing much more slowly re- 
cently, and hence one automatically finds that at high red- 
shift structure should be more advanced than in the critical- 
density case. Nevertheless, for completeness we perform here 
a calculation of the constraint arising from damped Ly- 
man alpha systems, following the strateg y based on Press- 
Schechter theory utilized by Liddle et al. (1996a) to extend 



the usual calculation to the open universe case. 

The contribution of the baryons present in damped Ly- 
man alpha systems to the mean cosmological mass density 
at some redshift z is given by (e.g. Padmanabhan 1993, p. 
347) 



Pdlas(z) = M m P {N) 



AN 
cdt 



(21) 



where is the baryonic mean molecular weight, m p is the 
proton mass, (N) is the mean H I column density as a func- 
tion of redshift and AN/ cdt is the number density of damped 
Lyman alpha systems per unit length interval at that red- 
shift. As dN/dz is a measurable quantity, the only depen- 
dence of pdlas(z) on the assumed cosmology comes from 
dz/cdt. Using the Friedmann equation it is straightforward 
to show that, for the cosmological models we are interested 
in, the amount of baryonic matter deduced from observation 
at a given z is given by the following function of S7o : 



^DLAS(^C), «) = fiDLAs(^o = 1, z) 



n 

Q(z) 



(22) 



For our cosmology, f2(z) is given by equation (J7J) , though 
the above result applies even to non-flat cosmologies with 
only matter plus a cosmological constant provided that the 
appropriate Q(z) is used. Expressing S7DLAs(no, z) in terms 
of its value assuming a critical-density cosmology is useful 
since most quoted observational results for Odlas^) assume 
such a cosmology. 

Estimates of the abundance of baryons in damped Ly- 
man alpha systems have recently b een p rovided at redshifts 
3 and 4 by Storrie-Lombardi et al. ( 19951) , the redshift 3 data 
giving a somewhat lower result than the m uch used earlier 
data of Lanzetta, Wolfe & Turnshek (1995). The constraints 
from the redshift 3 and 4 points are very similar, and we will 
use the latter. Using the above expression we then have 



go 

n(z = 4) ' 



(23) 



0dlas(2 = 4) = (0.0011 ± 0.0002) h' 1 

where the uncertainty is la. As mentioned before, Q(z) is 
given by equation (Q) and at redshift 4 is very close to unity 
for any reasonable f2o- Expression (^) conservatively as- 
sumes all the gas to be in the neutral state. 

Since we assume all the dark matter to be cold, a rea- 
sonable hypothesis is that the total density of the systems 
is larger than the baryonic density by a factor flo / Qb ■ Then 
the fraction of mass which is in bound objects of mass at 
least the typical mass M is given by 



/(> M, z = 4) > (0.069 ± 0.021) h 



n(z = 4) 



(24) 



where w e have interpreted the error on Q,b!i 2 in Copi et al. 
(19951) as 95 per cent confidence and added it in quadra- 
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ture. Using Press-Schechter theory (Press & Schechter 
197 4|), <| he theoretical prediction for this quantity is 



Frenk 1991; White et al. 1993). Typical quoted values sug- 



f(>M,z) 



Q(> M(R),z) 



erfc 



5c 



V2a{R, 



(25) 



where <5 C is a threshold parameter to be fixed via TV-body 
simulations and the smoothing to obtain a(R) is carried out 
via the top- hat window of equation (ji)) . 

The most conservative assumption is that the damped 
Lyman alpha systems have not fully collapsed along all their 
axes. This corresponds to a lower choice of thre shold than 
the usual 1.7. so. follow ing Liddle et al. (1996b), we adopt 
5 C = 1.5 (Monaco 1995). The minimum mass of the damped 
Lyman alpha systems (from the assumption that they even- 
tually give rise to rotationally su pported discs as indic ated 
by lower redshift observations (jLanzetta et al. 1995)) is 



taken as lO lo /i" 1 M flHaehnelt 1995| ) 



5 MISCELLANEOUS CONSTRAINTS 

In this section we consider other constraints which moti- 
vate the preferred choice of tto. They either fall outside the 
description of large-scale structure, or correspond to con- 
straints from large-scale structure that we have insufficient 
technology to compute ourselves and must take at face value 
from other papers. 



5.1 Age of the Universe 

The age of a flat universe is given by^ 



Hot = - 



3 VI - fin 



In 



i + yi - Qo 



(26) 



At fixed Ho, lowering flo increases the age. However, as al- 



ready noted for the case of an open universe ( Liddle et al 



1996a), observations such as the galaxy correlation function 
more or less fix the combination QoHo, rather than Ho it- 
self. With that combination fixed, lowering Qo actually leads 
to a reduction of the age of the model universe (though a 
less severe one than in the open case) . Consequently, ACDM 
models that fit other observations are younger the lower the 
density. 

Observational age determinations remain controversial, 
perhaps solely because the older estimates are in conflict 
with almost any cosmology (and most seriously with low- 
density cosmologies, as seen in our figur es) . However, a re- 
cent analysis by Chaboyer et al. (1996) (see also Bolte & 
Hogan 1995; Jiminez et al. 1996) suggests a 95 per cent 
lower limit of 12 Gyr, which is weak enough to leave most 
possibilities intact. 



5.2 Baryons in clusters 

The piece of evidence which is presently providing the 
strongest push for researchers to consider low-dens ity mod- 
els is the question of the baryon density in clusters ( White & 



[TjThere is an equivalent expression where the logarithm is re- 
placed by sinh -1 ^^(1 — Q.o)/£lo \ ■ 



gest that about 20 per cent of material in clusters is in the 
form of baryons, primarily as hot intracluster gas. Taken at 
face value, this is consistent with standard nucleosynthesis 
only if the total density is substantially less than critical. 

Although a pressing concern, there remain some po- 
tential problems with the interpretation of the observations 
(see for example the summary by Steigman & Felten 1995). 
A standard argument is that clusters are such large objects 
that the fractional baryon density within them must be rep- 
resentative of the Universe as a whole; however, individ- 
ual measurements of the fractional baryon density are not 
all consistent with each other on a cluster-by-cluster basis 
(e.g. White & Fabian 1995), so unless the observational er- 
rors have been somewhat underestimated this assumption 
must fail at some level. Nevertheless, this doesn't seem to 
be too important an uncertainty. For example, an unreason- 
ably crude method of bringing the White & Fabian (199E) 
results into agreement with one another is to throw out the 
clusters with the highest four or five determinations; even 
this only reduces the best-fitting baryon density by about 
10 per cent. Retaining all of their preferred sample of 13 
clusters, White & Fabian find 



Ob 
Oo 



0.05Z 



h 



-3/2 



(27) 



at the 95 per cent confidence level. A further question con- 
cerns whether or not the mass determinations (of either the 
gas or the total mass) are sufficiently sophist icated : for ex- 
ample, in a recent paper Gunn & Thomas ( 1996 ) carried 
out an analysis modelling the cluster gas as a multiphase 
medium, and found that the inferred fractional baryon den- 
sity can be reduced by a quarter or even a half. Finally, 
the status of standard nucleosynthesis theory has been less 
static recently than for many years, with the upper limit on 
the baryon density rising substantially. Copi et al. (1995a,b) 
advocate the central value Qsh 2 = 0.016 that we are adopt- 
ing, and their 95 per cent confidence upper limit is 0.024 
(see also Kernan & Krauss 1994; Sasselov & Goldwirth 1995; 
Krauss & Kernan 1995; Hata et al. 1995, and for an oppo- 
site opinion see Hogan 1995). Some recent deuterium abun- 
dance measurements in quasar absorption systems suggest 
even higher baryon f ractions (Tytler, Fan & Buries 1996 ; 
Buries fc Tytler 1996j). 



Ignoring mass determination un certai nties beyond 
those accounted for by White & Fabian (1995), in combina- 



tion the cluster and nucleosynthesis limits imply a central 
value f2o = 0.32/i~ 1/ ^ 2 . Taken at face value the error bars 
exclude the critical-density case for any reasonable h (the 
95 per cent upper limit being 0.53/i -1 / 2 ), but in our view , 
shared by Babul & Katz ( |l993[ ) and Gunn & Thomas ( ]l996| ) 
(see also Balland & Blanchard 1995), the theoretical uncer- 
tainties remain great enough that it would be extremely pre- 
mature to abandon critical-density models completely. It is 
clear however that present interpretation of this vital con- 
straint does favour values of Qo in the range between 0.2 
and 0.6. 



5.3 Velocity flows 

A direct analysis of the veloc ity flow data has been carrie d 
out by Nusser & Dekel ( 1993 ) and Bernardeau et al. ( 1995 ). 
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By analysing statistics of the velocity field in the weakly non- 
linear regime, they conclude that Qo must exceed 0.3 at the 
95 per cent confidence level. The same result has also been 
obtained from an ana lysis of velocity outflows in void regions 
( Dekel & Rees 1994 ) . Sum maries of these results are give n 
in reviews by Dekel (Il994|) and Strauss & Willick (h 995|) . 



The parameter region excluded by this constraint has only 
minimal overlap with our allowed regions, and always in a 
region where the age of the universe is anyway suspiciously 
low. So our analysis independently supports the conclusion 
of these papers. 

The results of those methods are quoted independently 
of any assumptions about the underlying power spectrum 
(though the former method does have some dependence on 
the assumption made). It would be expected that if the anal- 
yses were repeated with specific power spectra input, then 
the constraints would strengthen^]. 

5.4 Lensing limits on A 

An extremely comprehensive account of lensing limits on 
the cosmol ogical constant has recently been presented by 
Kochanek (1995). He concludes that the 95 per cent con- 
fidence limit for flat universes is Qo > 0.34, a constraint 
which is similar to that obtained from velocities. However, 
he notes that this limit is sensitive to several possible sys- 
tematic biases, and it can be relaxed to some extent by the 
introduction of extinction due to dust obscuration in the 
lensing population of mainly E/SO galaxies. 



5.5 Type la supernovae 

Prospects now look extremely good for applying a classic 
cosmological test, using type la supernovae as standard can- 
dles at very large distances in order to study the decelera- 
tion of the Universe. Preliminary results from the Supernova 
Cosmology Project (Perlmutter et al. 1996) favour a decel- 
erating Universe rather than an accelerating one, which if 
confirmed will give a stronger constraint than the lensing 
limits of the previous subsection. A decelerating flat universe 
requires Qo > 2/3. With many more supernovae already be- 
ing analysed, this is perhaps the most promising route of all 
towards constraining low-density flat universes. 



5.6 Galaxy counts 

Faint galaxy counts in principle offer a fairly direct probe of 
the geometry of the Universe to high redshifts, but in prac- 
tice the separation of geometry from evolution of the galaxy 
population ha s proven difficult. A recent analysis using HST 
observations (Driver et al. 1996) concludes with some defi- 
niteness that for a flat universe Qq must exceed 0.2, which is 
weaker than other constraints. A similar constraint had al- 
ready been obtained from deep number counts by Gardner, 

t After our paper was submitted, analyses based on computation 
of the power spectrum of mass fluctuations directly from the PO 



TENT Mark III cat alogue were produced ( [Kolatt fc Dekel 1995 
Zaroubi et al. 1996), including comparison with a range of cos 



mological models. Where they overlap, our results seem in good 
agreement with theirs. 



Cowie & Wainscoat ( 1993| ). However, Driver et al. (199f) 
add that with further assumptions concerning correctness 
of current morphological classifications, that limit can al- 
ready be dramatically strengthened to Qo > 0.8 for a flat 
geometry. Further development of this constraint is clearly 
a promising avenue. 



6 RESULTS 

Our results are shown in Figs 1 and 2, which indicate the 
microwave background and large-scale structure constraints, 
along with some age contours for reference. Let us first dis- 
cuss Fig. 1, which shows the situation without gravitational 
waves. 

Fig. 1 shows that incorporating the possibility of tilt 
greatly increases the space of viable parameters. By far the 
most important constraint is the shape of the galaxy corre- 
lation function; this constraint is reinforced by the cluster 
abundance constraint which lies in a similar region but is 
typically weaker. For low n, the allowed region is then fur- 
ther trimmed by the damped Lyman alpha system abun- 
dance. The POTENT constraint is weaker than the others 
in all parts of parameter space and we omit it for clarity. At 
low Qo, models start to become dangerously young. 

One can go to very high n in these models. Although 
we haven't indicated it in a plot, we have found that there 
remains an allowed region even at n — 1.4. However, once 
n > 1.3, a strong anti-bias is necessary unless h is well below 
0.5. Also, such models predict a huge Doppler peak, and 
require a redshift of reionization which is just right to reduce 
the large peak down to the standard size. 

In almost all regions of parameter space the height of 
the Doppler peak is not constraining. The exception is when 
n — 0.8, where it begins, quite dramatically, to cut off the 
high-h, region. To avoid confusion through plotting yet an- 
other line, we have indicated this simply by not shading the 
region (centred around Qo = 0.5, h — 0.75) that is allowed 
by other constraints but not the Doppler peak height. The 
exclusion is because for these parameters the peak is not 
high enough to fit the observations. As n is decreased, the 
Doppler peak constraint rapidly becomes very strong; al- 
ready by n = 0.75 there is no allowed region at all. Note 
though that the value of n at which this happens would be 
changed were the baryon density allowed to be substantially 
higher than indicated by nucleosynthesis. 

The main problem with the Doppler peak constraint 
is that the upper limit cannot be used due to the possi- 
bility that reionization might bring the peak back down to 
an acceptable height. However, some regions of otherwise 
permitted parameters do require specific assumptions about 
reionization and we will discuss this shortly. 

Except in rather extreme circumstances, the region al- 
lowed by these large-scale structure constraints is not fur- 
ther reduced by including the miscellaneous constraints from 
bulk velocity flows, gravitational lensing and galaxy counts 
that we discussed in Section 5, which favour Qq > 0.3. 
Rather, our results provide independent support to the con- 
clusions of those analyses. Type la supernovae ( perlmutter 
et al. 1996) seem the most promising route to a stronger 
constraint in the immediate future. 

In Fig. 1, we have shaded the allowed regions in two 
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Figure 1. The constraints plotted in the f2o— h plane, for different n in the case without gravitational waves. All constraints are plotted 
at 95 per cent confidence. The models are normalized to the COBE data. The solid lines are limits from the shape parameter, the 
dashed lines from the cluster abundance and the dot-dashed lines from the abundance of damped Lyman alpha systems. The POTENT 
constraint is not shown as it is weaker than all the other constraints. Contours of constant age are shown as dotted lines. The allowed 
region is shown with two different shadings, both highlighting the parameter space not excluded at more than 95 per cent confidence on 
any single piece of data. The lighter shading shows models where the optical galaxies have to be anti-biased at 8 h -1 Mpc. The unshaded 
region in the n = 0.8 plot which is allowed by all plotted data is excluded by Doppler peak height. 



different styles, to indicate whether anti-biasing of optical 
galaxies at 8/i _1 Mpc is required or not. In doing this, 
we take the COBE normalization to its 95 per cent lower 
limit, so the light shaded regions definitely require anti- 
biasing. Whether this is excluded on any particular observa- 
tional grounds is unclear. Motivated by Klypin et al. (1995), 
we have also tried to calculate to what extent anti-bias is 
required at scales presently in the non- linear regime. We 
choose to consider the APM power spectrum measurement 
(Baugh & Efstathiou 1994) at lh^ 1 Mpc. Using analytic 



approximations to estimate the non-linear power spectrum 
( iHamilton et al. 199l|; Peacock k Dodds 1994 pain, Mc 
fc Whi |e 19951 [Peacock fc Dodds 1996| ) for the models we 
are interested in, we can determine which require an anti- 
bias at this scale. As it turns out, the resulting constraint 
is not much stronger than assuming no anti-bias at Sh' 1 
Mpc using the linear power spectrum. In any case it seems 
less clear that anti-biasing is physically forbidden so far into 
the non-linear regime. The question of whether anti-bias is 
permitted will perhaps be accessible to new hydrodynam- 
ical A^-body simulations, which might indicate whether or 



not forbidding anti-bias at 8 h 1 Mpc can really be adopted 
as a firm constraint^. 

Our conclusion in the case of no gravitational waves is 
that the large-scale structure data alone permit Qo as low as 
about 0.30, and show no preference for any particular value 
above that. However, at critical density the required Hubble 
parameter is very low ( though see the more extensive dis- 
cussion in Liddle et al. (1996b) which also considers adding 
a hot dark matter component). The ACDM model is com- 
patible with observations for a wide range of values of the 
tilt. 

Fig. 2 shows what happens if gravitational waves are in- 
troduced in accordance with the power-law inflation model. 
This model requires n < 1 so we only show two values. The 
introduction of gravitational waves leads to a much smaller 
allowed region than in their absence, with n = 0.9 the small- 
est allowed value. At n — 0.9, the shape and cluster con- 
straints leave only a very narrow strip, which is then totally 
excluded by the Doppler peak constraint. For n — 0.95 there 



S Also see Kauffmann, Nusser &; Steinmetz ( L99J ) for a recent 
discussion of the physical origin of biasing; their analysis does 
not favour anti-biasing. 
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h 0-6 




h 0.6 




Figure 2. As Fig. 1, but with gravitational waves as given by power-law inflation. In each case where only a single line appears, the 
allowed region is to the upper right. In both of these plots, the unshaded region which is allowed by all plotted data is excluded by 
Doppler peak height. 



is an allowed region. The limit n > 0.9 in the power-law 
i nflatio n model is stronger than that given by Liddle et al. 

in the critical-density case, because here we consider 



( [!996bD 



the additional constraint coming from the Doppler peak. 

Although we will not fully investigate the issue here, 
the Doppler peak constraint does have some implications 
for reionization. Because ACDM models at low density have 
a very high normalization (especially at high redshift), early 
reionizatio n is more likely in th ese models than in other cos- 
mologies ( Liddle fc Lyth 1995 ). However, the physical un- 
derstanding of the reionization process remains poor enough 
that most estimates of the reionization redshift can be in- 
terpreted only as upper limits. 

We'll restrict ourselves to the case with no gravitational 
waves, and for simplicity assume the central nucleosynthe- 
sis value for Qb (to which the results are quite sensitive). 
In some regions of parameter space, any early reionization 
would remove the Doppler peak sufficiently to be incom- 
patible with observation, while in others early reionization 
is necessary to bring it down to acceptable levels. Using 
the formulae given earlier, we find the following results. For 
n > 1.1, reionization at any 2r less than 15 is not enough 
to allow any interesting models. For lower n, reionization 
becomes less necessary in order to lower the Doppler peak; 
indeed, it becomes a worry that it may lower it too much. If 
n = 0.8, then zr > 15 is sufficient to exclude all the other- 
wise allowed region. A full examination of the relationship 
between large-scale structure constraints and the required 
reionization properties, taking into account the uncertainty 
in the baryon density, is clearly warranted. 



7 CONCLUSIONS 

Our results extend previous work both in terms of the range 
of parameter space studied and in terms of the observa- 
tions considered. By including the possibility of tilt and 
gravitational waves, we have explored the parameter space 
available for inflationary-based ACDM models. Recent work 
(Krauss & Turner 1995; Ostriker & Steinhardt 1995; Bagla 
et al. 1995) has favoured flat low-density models primarily 



not for large-scale structure reasons, but rather as a resolu- 
tion of two separate conflicts, the first being age verses the 
Hubble constant and the second being the cluster baryon 
fraction verses nucleosynthesis. We have shown that while 
compatibility with large-scale structure does not in itself se- 
lect favoured values of flo, it is compatible with low-density 
ACDM models. 

It is interesting to compare these results with the open 
CD M case explored using the same techniques in Liddle et 
al. (1996a). The principal difference is that at low fio the 
COBE normalization is very different, so that in the open 
case anti-bias is not necessary. However, other observations, 
such as the cluster abundance, shift as well to maintain a 
substantial allowed region. The allowed regions are fairly 
similar in the ACDM and open CDM cases, with no partic- 
ular preference for either on present data, though the open 
models tend to prefer a slightly higher value of n. In going 
to a ACDM model, there is only a very modest gain in age 
relative to an open CDM model, so the principal motivation 
for introducing the cosmological constant is in maintaining 
consistency with standard inflationary models rather than 
because of an age crisis. 

Returning to our discussion of the ACDM model, an 
extremely useful result is the fitting function, equation (hot) , 
which allows a quick normalization of all A models, since to 
an excellent approximation it is independent of the parame- 
ters h and Qb and the presence of any component of hot dark 
matter. We have also made the first detailed comparison 
of these models with the observational situation concerning 
intermediate-scale CMB anisotropies; though observations 
there are still at a very primitive level compared with what 
one hopes will become possible in the future, already one can 
obtain new constraints against regions of parameter space 
involving substantial gravitational waves, which are viable 
against all other present observations. We have also noted 
that certain areas of parameter space are viable against this 
constraint only if reionization occurred suitably; in some re- 
gions it is required to happen fairly early, while in others it 
is forbidden from doing so. 

In the case of a scale-invariant spectrum with no gravi- 
tational waves, our results are in broad agreement with those 
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obtained by other authors. The main challenge from large- 
scale structure for the low-density versions of these models 
continues to be the question of whether it is reasonable to 
believe that optical galaxies are un- or even anti-biased. We 
have included more constraints than the earlier treatments 
but for scale-invariant spectra they are not more constrain- 
ing. 

By extending the parameter space to include tilt and 
gravitational waves, we have found a large region compati- 
ble with present data. Amongst low-density models, it seems 
Qo in the range 0.4 to 0.5 is better than Qo = 0.3 (continu- 
ing a historical upward trend in the preferred ' low-density' 
Qo), as is also the case for open CDM models (Liddle et al 
1996aQ 



particu 



However, large-scale structure studies do not look 



irly promiaing for lifting the degeneracy of allowed 



Qo values; the types of evidence we discussed in Section 5 



appear prore promising for delivering ail 
or agai nst the coamological constant. 



ultimate verdict for 
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